CP Violation in neutrino oscillation and leptogenesis 
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We study the correlation between CP violation in neutrino 
oscillations and leptogenesis in the framework with two heavy 
Majorana neutrinos and three light neutrinos. Among three 
unremovable CP phases, a heavy Majorana phase contributes 
to leptogenesis. We show how the heavy Majorana phase 
contributes to Jarlskog determinant J as well as neutrinoless 
double P decay by identifying a low energy CP violating phase 
which signals the CP violating phase for leptogenesis. For 
some specific cases of the Dirac mass term of neutrinos, a 
direct relation between lepton number asymmetry and J is 
obtained. For the most general case of the framework, we 
study the effect on J coming from the phases which are not 
related to leptogenesis, and also show how the correlation can 
be lost in the presence of those phases. 
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Finding any relation between baryogenesis via lepto- 
genesis [1] and low energy CP violation observed in the 
laboratory is a very interesting issue [2]. The CP viola- 
tion required for leptogenesis stems from the CP phases 
in the heavy Majorana sector, whereas CP violation mea- 
surable from the neutrino oscillations [3] can be described 
by the neutrino mixing matrix. One interesting question 
concerned with the low energy leptonic CP violation is 
whether it can be affected by the CP violating phases re- 
sponsible for leptogenesis. Several people [4] have already 
discussed some potential connections between low energy 
CP violation and leptogenesis by using some ansatz, but 
it is still unclear how large the former can affect the lat- 
ter in general. The major difficulty to quantify such a 
connection occurs due to lack of the available low energy 
data to fix parameters of the seesaw model. 

The purpose of this paper is to examine in a rather gen- 
eral framework how leptogenesis can be related to the low 
energy CP violation by determining the parameters as 
many as possible from available low energy experimental 
results and cosmological observations. In order to make 
a quantitative analysis of the connection between low en- 
ergy leptonic CP violation and leptogenesis, we consider 



the minimal CP violating seesaw model which has two 
heavy Majorana neutrinos and three light left-handed 
neutrinos; (3,2) seesaw model. As will be shown later, to 
break CP symmetry, the required minimal number of sin- 
glet heavy Majorana neutrino is two in the seesaw model 
with three light lepton doublets. This (3,2) seesaw model 
is consistent with recent data of neutrino oscillations and 
contains 8 real parameters and 3 CP violating phases in 
the neutrino sectors which make this model more con- 
strained and predictive compared with the general (3,3) 
seesaw model [5] with 18 parameters. We will show that 
while all three CP violating phases contribute to low en- 
ergy leptonic CP violation, only a single CP violating 
phase contributes to leptogenesis. We will also investi- 
gate how large the CP phase responsible for leptogenesis 
contributes to low energy CP violation by determining 
the independent parameters from available experimental 
results and cosmological observations. Finally, we will 
discuss the potential implication of CP violation measur- 
able from neutrino oscillations on leptogenesis. 

Let us begin our study by considering the leptonic sec- 
tor of the (3,2) seesaw model. In a basis where both 
heavy Majorana and charged lepton mass matrices are 
real diagonal, the Lagrangian is given by: 

1- 



C = -liLmiikR - VLiruDijNRj - -{NrjYMjNrj, (1) 

where i — 1,2, 3, j — 1,2 and the Dirac mass term m/^ 
is 3 X 2 matrix. Here, we remark that the Dirac mass 
matrix mu contains 3N — 3 unremovable CP phases if 
we take N singlet heavy Majorana neutrinos in this basis. 
Thus, one can easily see that at least two singlet heavy 
Majorana neutrinos are required to break CP symmetry 
in the seesaw model with three lepton SU{2) doublets. 
The 3x2 matrix mu can be generally parameterized as: 



mo = Ul 




(2) 



with, 
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Ul = O23(^L23)C/i3(0L13,<5l)Oi2(^L12)P(7l), 

/ PriQ Hn Qin Hr> \ I '' 

Vr = 



cos 9r sin 9r 
— sin 6r cos 9r 











(3) 



where Oy and Uij 



P{1l) 



denote the rotations of (i,j) plane, 
diag.[l, e"*'''^/^, e'^^/^], and m,2 and TO3 are real 



and positive. Without loss of generality, we can choose 
< ni3. The allowed range for the angles and the 
phases is [— 7r,7r]. There are three CP violating phases, 
which appears in Vr, and Sj^ in Ul- In a different 
basis with complex Mj, can be interpreted as a heavy 
Majorana phase. The lepton number asymmetry for the 
lightest heavy Majorana neutrino (A''i) decays into 
[7] is given by; 



Ti-ri 

£l = 



3Mi I m[{{m 0^1710)12} 



(4) 



where V = \/Attv with v = 246 GeV, and 



Im[{{mD^mD)i2y 



sm'^20Rsm2-fR. (5) 



We see that CP violation concerned with leptogenesis 
can be possible only if the mixing angle Or and CP vi- 
olating phases -jr for the heavy Majorana neutrinos are 
non-zero. For our purpose, let us now study how the 
phase 7fl contributes to CP violation in the neutrino os- 
cillations which is usually described in terms of the MNS 
neutrino mixing matrix [6]. The effective mass matrix 
for hght neutrinos is given by rueff = —mD-^fno^ = 
—ULrnVR^VR^m^U]^, and is diagonalized by the MNS 
mixing matrix as U\j^gm,effU1jjqg = diag[ni,n2,n^]. 
Then, the MNS mixing matrix is decomposed into two 
mixing matrices as follows; 



Umns = UlKr, 



(6) 



where Kr is a unitary matrix diagonalizing the matrix 
Z = —rnVR^VRrn^ and parameterized by, 

1 \ / 1 

Kr = \ COS0 sinOe-^'f' e*" ) . (7) 
-sin6'e** cos6' / \0 e"* 

Then, -Kl^rnVR^V^rn^ K'^ = diap.fO, n2, ns]. Note 
that the (3,2) seesaw model predicts one massless neu- 
trino. In addition, 9, 4> and a are determined as: 



tan 26* = 



Arg.{Z22Z2i, + Z23Z33), 

21-^22-^23 + -^23-^33! 



(8) 



|^33P - \Z22\^ ' 

2a = Arg.[cos^ eZ22 + sin^ 61^336"^'"^ - sin26'Z23e-''^]. 

We remark that the mixing angle 9r and CP violating 
phase 7_R have been transferred io9,(j) and a. As one can 
see from the above formulae, leptogenesis occurs only 
if the mixing angle 9r and CP violating phase ^r arc 
non-zero, which in turn implies non-vanishing (j), a, 
via Eq.(8). As we will see below, the CP phase (j) con- 
tributes to CP violation in neutrino oscillations, so that 
it is anticipated that there is correlation between CP vio- 
lation generated from neutrino mixings and leptogenesis. 



To see this concretely, let us compute Jarlskog deter- 
minant [8] J = Im[UMNSelUMNSe2UMNSfilUMNSfi2], 

which is proportional to the CP asymmetry in neu- 



trino oscillation, AP 
4J (sin(^ 



£) + sin(^^)+sin(^ 



#^))- By 



usmg 



2E 

Eq.(6), we obtain: 

J = sin 26li2 sin 2^1,13 [cli3 cos 26 sin 5l sin 2^i23 

o 

+ CL12 sm26sm{5L - - <l>) cos 2^1,23 

- ^31,1251,13 sin 2^ sin 2^1,23 sin(25i:, - Jl - (t>)] 

+ ^ sin 26 sin 2^1,23 sin(7i + (p) x 
o 

{sm26Li2CLi3SLi2 - sm 26 L13S L13CL12) ■ (9) 

Prom the expression of J, it is obvious that all three CP 

violating phases , 7l and (f> contribute to CP violation 
in the neutrino oscillations, and that the CP phase 
always hangs around (f). Since only (f) is closely related to 
leptogenesis, in order to investigate the interplay between 
CP violation for leptogenesis and low energy leptonic CP 
violation, we should determine the contributions of the 
phases ((5l,7l) and (j) separately as well as to fix the 
parameters 9's. 

Before discussing the correlation between both CP vi- 
olations, let us study how we can get some information 
on the mixing angles and CP phases from the available 
experimental and cosmological results. The mixing an- 
gles and CP phases can be classified into two categories, 
one contains 9, (j) and a which are related to phenomena 
at high energy and the other contains parameters in Ul- 
First of all, we show how wc can estimate the allowed 
values of CP violating phase </> and mixing angle 0. The 
information on </> and a may come from the constraints on 
light neutrino mass spectra as well as cosmological con- 
dition for leptogenesis. To see this, we first present the 
parameters 7fl,^fl,rn2,m3 and lepton number asymme- 
try El in terms of some physical quantities which will be 
taken as inputs in numerical calculation. Here, we choose 
the heavy Majorana neutrino masses (Mi, M2), their de- 
cay widths (Fi, r2), and light neutrino masses (n2, 713) as 
the physical input parameters. As will be clear later, it 
is convenient to define two parameters a;i(i = 1, 2); 

(mD^mD)j 



Xi 



M, 



(10) 



Then, by considering the light neutrino mass eigenvalue 
equation, det[m(,ffrn^^^^ — n^] = 0, the lepton number 
asymmetry ei and the phase 7^ can be written in terms 
of x\,X2,n2, and 713, 

3Mi 



£1 = - 



4a;il/2 



^{{n.f-{x.f){{x+f-{n+f), (11) 



cos 2'yR 



Ttri X-t 3jO 



2{xiX2 — n2n3) 



(12) 



2 



where n± = ns ± 712 and x± = xi ± X2. There are 
two solutions of Eq.(12) leading to negative ei; and 
7ii — TT for < 7fl < n/2, which in turn gives posi- 
tive baryon number via sphaleron process. Next, let us 
present the parameters, 9r, 7712 and ma in terms of the 
above 6 physical quantities. Prom the eigenvalue equa- 
tions for Vfi^^^^==-Vji^ we can express 6ji, m^, and ms 
as follows; 



{ml,m\ 



\)=^/W[M2{cj+-p,cj++p), (13) 
(cos^H,sin^«) = {,[^.-,[^^^ , (14) 



2p 



where a± = P = {xiX2 - nin^,) + a'^_, and 

R = M1/M2. We also determine (f>, a, and 9 with a given 
set of parameters (xi, 3:2, i?, 712, ria) by using the same 
procedure given in Eqs.(7,8). We take R = 0.1. In order 
to determine the values of cj), 9 and a, it is necessary to 
determine those of Xi{i = 1, 2). Let us now show how the 
variables Xi can be constrained. Prom the neutrino mass 
eigenvalue equation, it follows that 



Fi - X2\ < na - n2, 
na + n2 < xi + X2- 



(15) 



Prom the experimental results for the neutrino oscilla- 
tion, let us take = \J Am^^^ ^ 5 x 10-2 eV and 
"2 = V^m'ioiar = 7 X 10-3 eV (LMA) [9]. FromEq.(15), 
the lower bound on xi is 0.007 eV. By solving the Boltz- 
mann equation [11], we can obtain a value of Yl = — , 
i.e., lepton number density {ul) normalized by entropy 
density (s). When solving the Boltzmann equation, we 
need the value of Si. For fixed xi, one can get the max- 
imum value of —Si which gives the maximum — Yl via 
Boltzmann equation. In Pig.l, we plot the maximum 
lepton number density — Yl predicted from Eq.(ll) as 
a function of xi for several fixed Mi .We set the initial 
conditions for Boltzmann equation at 10^^ GeV and we 
take the distribution of the heavy majorana particle Ni 
in thermal equilibrium and Yl = at the temperature. 




FIG. 1. The maximum possible lepton number density 
(— Yl) as a function of Xi for three different Mi. From up 
to bottom, Ml = 3, 2, 1 X lO" GeV. 



The allowed values of — Y^ consistent with baryogen- 
esis are presented by shaded band in Fig.l. Thus, we 
can obtain the allowed region of xi for a fixed Mi . How- 
ever, there is no allowed value of xi for a rather lower 
value of Ml < 1.0 x 10^^ GeV, which in turn leads to the 
lower bound on Mi. By using the allowed region for xi 
as given in the above, we can estimate the allowed region 
of X2 via Eq.(ll) again. Figure 2 shows how we can get 
the allowed region of X2- Por example, for a given set 
Ml = 2 X 10" GeV, and xi = 0.03 eV, we obtain the 
allowed range for X2 as 0.03 < X2 < 0.07 eV. 
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FIG. 2. Lepton number density (— Yl) as a function of X2 
for Ml = 2 X 10" GeV. Solid line corresponds to xi = 0.01 
eV. The long dashed line corresponds to xi = 0.03 eV and 
the dotted line corresponds to xi = 0.05 eV. 



Let us move to the other category of the parameters, 
9li2,9l23,9ls13,1l,Sl in Ul, which are not determined 
from high energy phenomena, but must be related to the 
low energy MNS mixing matrix, Thus two of them can 
be determined from the neutrino oscillation experimental 
results. For simplicity, we focus on the case with the 
small mixing angles and 9, which is consistent with 
Chooz experiment [10]. In this case, the MNS mixing 
matrix is given, in the leading order, by 



Umns — 

CL12 SL12 

— SL12CL23 CL12CL23 
SL12SL23 —CL12SL23 

xP{a',-a'). 



' + SLi2See 

SL23 
CL23 



-iff}' 



(16) 



where (j)' = (j) + 7^, a' = a ^ ^ and P{a' , —a') = 
diag.[l, e'" , e~'^" ]. Note that we do not present the sub- 
leading contributions in {UMNs)ij, {ij) 7^ (e3), which are 
comparable to {UMNs)e3- Taking f?23 — 9i2 ~ ^ which 
lead to bi-large mixing pattern, in this approximation, 
{UMNs)e3-, J and |(me//)ee| are given by: 



|(t^MAfS)e3 



SLiae 



-i<5 



\/2 



J ~ ^ (sli3 sin 5l + ^ sin ( 
4 V V2 



{meff)ee\^\^e'"^' 



30 



+ n3(5Li3e-*^- + -|e-^^')2|. (17) 



In principle, we are able to fix three unknown parame- 
ters; 7L,Sii3 and 5l once the left-hand sides of Eq.(17) 



3 



are measured. It is then possible to quantitatively see 
whether the low energy CP violation denoted by J is 
dominated by leptogenesis phase </> or by the CP vio- 
lating phases 7^ and Sl which are not related to lep- 
togenesis. First of all, let us study the interesting case 
of 9li3 = 0, which makes the analysis more predictive 
because a CP violating phase Sl simultaneously sup- 
pressed. This can be understood as the extreme case of 
sin^iia <C sin 6*. Interestingly enough, this case dictates 
that the origin of |(C/j\/Ars)e3| niay come from the mix- 
ing angle 9 which is related to heavy Majorana neutrino 
sector. Jarlskog determinant is then simply given by. 



J = sin 20 sine 



1 



8V2' 



(18) 



Only 7i is a completely arbitrary parameter in this case 
and thus we can easily investigate how J can be affected 
by 71,. In other word, 7^ can be estimated through J 
in this case. If 7^ is turned out to be much smaller 
than (j), the measurement of CP violation in low energy 
experiment may directly indicate leptogenesis. 




epl * 10"6 

FIG. 3. Correlation between lepton number asymmetry ei 
and neutrino CP violation J, for Mi = 2 x lO" GeV. The 
large (small) contour corresponds to xi = 0.03(0.01) eV. The 
thick solid lines show the allowed region from baryogenesis. 



In Fig. 3, we show the correlation between lepton num- 
ber asymmetry ei and J in Eq.(18) with 7^ = 0. In 
each contour. Mi and xi are fixed and X2 is varied. For 
xi = 0.01 eV, we obtain \{UMNs)e3\ = ^ ~ 0.066 
and \{meff)ee\ = 0.008 eV, while for xi = 0.03 eV, 
|(C^AfJVs)e3| is in the range [0.15,0.2] and \{meff)ee\ is 
in the range [0.014, 0.018] eV. On the other hand, the 
case with sin 6 <^ sin 6'l 13, we see from Eq.(17) that J 
mainly depends on Sl , which has nothing to do with lep- 
togenesis. 

We have studied how CP violation responsible for 
baryogenesis manifests itself in MNS matrix and Jarl- 
skog determinant which signals low energy CP violation 
in neutrino oscillation. We have obtained cosmological 
constraints on CP violation and mixing which originate 
from high energy phenomena. Then using the low energy 
constraints we have showed it is possible to estimate the 
size and sign of baryon number in the most general case 
once \ {UMNs)e3\, neutrinoless double beta decay and CP 



violation of neutrino oscillation are measured. In a spe- 
cific case of the framework, a correlation between CP vi- 
olation in neutrino oscillation and leptogenesis has been 
studied and the size of J has been estimated. 
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